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Unpublished investigations in this laboratory have shown that Clostrid- 
tum kluyvert can metabolize acetic acid and ethanol under anaerobic con- 
ditions producing butyric and caproic acids in accordance with equations. 
(1) and (2). 


CH;CH,0OH + CH;COOH = CH;CH:CH:COOH + H;0 (1) 
2CH;CH20OH + CH;COOH = CH;(CH2)4COOH + 2H20 = (2) 


The relative yields of the two acids are determined by the relative amounts 
of acetic acid and ethanol available. If acetic acid is present in excess, a 
considerable amount of butyric acid is formed, while if ethanol is in excess, 
caproic acid is the main product. These relations suggest that butyric 
acid may be an intermediate in the synthesis of caproic acid from acetic 
acid. This is supported by the observation that ethanol and butyric acid 
can be converted to caproic acid according to equation (3). 


CH;CH,0H + CH;CH:CH,COOH = CH;(CH:),COOH + H:0 (3) 


In the present investigation we have obtained positive proof for the con- 
version of acetic acid to caproic acid via butyric acid by studying the action 
of Cl. kluyvert on media containing synthetic fatty acids labeled with the 
long-lived radioactive carbon isotope C'*. In addition we have obtained 
various types of evidence which help to elucidate the mechanism of fatty 
acid synthesis from C, molecules. 

Experimental Methods.—Carboxyl-labeled acetic, butyric and caproic 
acids were prepared from C*O, and methyl iodide, propyl bromide and 
amyl bromide, respectively, by the Grignard synthesis.? 

The general technique of estimating the C'4 content of carbon dioxide 
and fatty acids has already been described elsewhere.? The barium salts 
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of the fatty acids dried at 100°C. were used in the radio-assay. The indi- 
vidual acids were separated from mixtures by the distillation procedure of 
Schicktanz, et al.* The purity of each acid was established by Duclaux 
distillation. The methods used in locating C'* in the individual carbon 
atoms of acetic and butyric acids have been described previously.” * ® 
The occurrence of C'4 in the carboxyl group of caproic acid was deter- 
mined by decarboxylating the barium salt. Control experiments with syn- 
thetic caproic acid showed that this method yields 52 per cent of the car- 
boxyl carbon as barium carbonate. 

The bacteria were grown in media of the following composition in g. per 
100 ml.: ethanol 0.3-0.8, sodium salt of labeled acetic or butyric acid 
0.3-0.6, yeast autolysate 0.3, pH 7.2, phosphate buffer 17/40, (NH,4)2SO, 
0.1, MgSO«-7H2O 0.01, FeSO..7H:O 0.0001, sodium thioglycollate 0.05, 


TABLE 1 
THE FERMENTATION OF CH;CH:OH anp CH;C*OOH sy C1. kluyveri 


Experiment 1 


COMPOUND MM cts./M1In./MM TOTAL CTS./MIN. 

‘Before growth 

Ethanol ca. 1.1 acs ae 

Acetic acid 0.830 18,900 = 300 15,700 

Carbon dioxide ca. 0.2 pata Sis. 
After growth 

Acetic acid 0.395 7330 + 150 2895 

Butyric acid 0.196 18,200 + 300 3570 

Caproic acid 0.257 24,500 + 500 6300 

Carbon dioxide 0.185 26 + 50 5 

Ethanol 0.217 260 = 50 74 

Non-volatile cpds. Ma ae 193 
Percentage recovery of C¥4........... pial s ciwte bis bale a nloteteleed be Bioh ave wea. 83 


traces of salts of Ca, Mn and Mo, distilled water. In some experiments 
sodium carbonate (0.025 g. per 100 ml.) was added as a sterile solution 
after autoclaving the medium. Oxygen was removed by means of either a 
pyrogallol-potassium carbonate or an Oxsorbent seal. The cultures were 
incubated at 32-35° until growth ceased before being analyzed. 

Results —When Cl. kluyveri was allowed to grow in a medium containing 
approximately equivalent amounts of ethanol and synthetic carboxyl- 
labeled sodium acetate, radioactive butyric and caproic acids were formed 
(table 1). The isotope was proved to be present in these acids, rather than 
in some associated compound, by establishing the constancy of the specific 
activities of the barium salts prepared from different fractions of a Duclaux 
distillation (table 2). A little C'* was also found in ethanol and unidenti- 
fied non-volatile compounds, but no significant amount was present in 
carbon dioxide. It should be noted that there is no net production of 
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carbon dioxide in this fermentation. The recovery of C' in all forms after 
growth was about 83 per cent of that initially added as acetic acid. 

A particularly significant fact revealed by this experiment is the de- 
crease in the molar activity of acetic acid during the fermentation. At the 
beginning of the experiment the activity expressed in cts./min./mM was 
18.9 X 10* while at the end the value was 7.3 X 10%. The simplest and 
most likely explanation for this reduction in activity is the oxidation of 
ethanol to inactive acetic acid or some closely related compound in equi- 
librium with it. This would cause a dilution of the labeled acetic acid. 


TABLE 2 
SPECIFIC ACTIVITIES OF BARIUM SALTS OF BUTYRIC AND Caproic ACIDS PREPARED FROM 
DucLaux DISTILLATIONS 
Experiment 1 (total volume = 110 ml.) 
FRACTION cTs./MIN./MG. OF Ba SALT 


Butyric acid 


0-20 ml. 47.5 
20-40 ml. 47.6 
40-110 ml. 46.4 

Caproic acid 

0-20 ml. 63.7 

20-110 ml. 62.2 


Another possible explanation for the decrease in activity of acetic acid 
is that it is in chemical and isotopic equilibrium with the higher fatty acids 
formed from it and ethanol. This was tested by two experiments in one 
of which butyric acid labeled in the carboxyl group was added at the be- 
ginning of the fermentation; in the second experiment butyric acid labeled 
on all four carbon atoms was used. In the experiment with carboxyl- 
labeled butyric acid (table 4), no C' could be detected in the final acetic 
acid, thus indicating that no exchange occurred between the carboxyl 
groups of butyric acid and acetic acid. In the second experiment with com- 
pletely labeled butyric acid, a little C'* was recovered in the final acetic 
acid but it was less than 9 per cent of the amount to be expected if butyric 
and acetic acids were in isotopic equilibrium, and even this small effect 
may have been due to a slight contamination of the original butyric acid 
by labeled acetic acid. An exchange of C!* between acetic acid and the 
higher fatty acids can therefore be excluded as an explanation for the 
observed decrease in activity of acetic acid. 

A third possible explanation for this effect is an exchange of C'* between 
ethanol and acetic acid as a result of reactions involving acetaldehyde or a 
similar compound as an intermediate. If such a reaction occurs, the 
residual ethanol should have the same molar activity as the acetic acid. 
Actually, the activity of the alcohol is of a lower order of magnitude 
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(table 1). Therefore an exchange of C'* between acetic acid and ethanol is 
not consistent with the observed results. 

Whatever the actual mechanism of the dilution of the active acetic acid, 
the magnitude of the effect is about what would be expected if all the 
ethanol is oxidized to acetic acid while the latter is being converted to 
higher fatty acids. To calculate the expected decrease of acetic acid 
activity on this basis, let us assume that (1) the fermentation proceeds in 
two successive steps represented by equations (1) and (3), and that during 
_ the fermentation all the ethanol is first converted to acetic acid or some 
compound in isotopic equilibrium with it. In step 1, the formation of 
butyric acid from ethanol and acetic acid, let x = the amount of C™ present 
in acetic acid at any time, x) = the initial C'*, x, = the C' at the end of 
step 1, V = the amount of acetic acid present at any time, Vo = the initial 
acetic acid, and V;, = the final acetic acid. Now during step 1, when a 
small amount (AV) of alcohol is oxidized to inactive acetic acid, the re- 
moval of active acetic acid by conversion to butyric acid will be equal to 
2AV. The loss of C' from acetic acid is therefore 


— Ax = —2AV/V — AV-x. (4) 
In the limit as AV approaches zero 
dx/x = 2dV/V. (4a) 
Integrating between the limits of Vp and V;, and x» and x; and changing to 
logio we get 
log x:/x%9 = 2 log V;/Vo (5) 
or 
x1 = x0(V;/Vo)?. (5a) 


In step 2 (equation 3), the conversion of alcohol and butyric acid to caproic 
acid, there is no change in the total quantity of acetic acid; V; is therefore 
aconstant. CC!‘ will nevertheless be lost from acetic acid if, as is assumed, 
it is an intermediate in the utilization of ethanol. When a small amount of 
alcohol (AA) is used, the loss of C'4 from acetic acid is 


— Ax = AA/(V, + MA)-x. (6) 
In the limit 
—dx/x = dA/V;. (6a) 


Integrating between the limits x, and x,;, and 0 and A,, where A, is the total 
amount of alcohol used in step 2 


In X1/X;z = 2.3 log x1/Xyz = Aj/V; (7) 


or 
log x, = log x1 — A;/2.3V, (7a) 





Oo VE = tee OO D 


Ye —_— > cP 
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When this method of calculation is applied to the data of table 1 where 
x0 = 15,600 cts./min., Vo = 0.830 mM, V; = 0.395 mM and A; = 0.257 
mM, it is found that x, = 3540 cts./min., and x, = 1840 cts./min. The 
calculated molar activity of acetic acid at the end of the fermentation is 
therefore 1840 cts./min./0.395 mM = 4660 cts./min./mM. This value is 
of the same order of magnitude though somewhat smaller than the observed 
molar activity of 7330 cts./min./mM. The discrepancy may be due to the 
fact that butyric and caproic acids are formed simultaneously rather than 
successively during part of the fermentation; this would raise the calculated 
activity of the final acetic acid and bring it into closer agreement with the 
observed value. The general conclusion to be drawn from a comparison of 
the calculated and observed decrease in the molar activity of acetic acid is 
that the results support the view that alcohol is oxidized to acetic acid or a 
closely related compound which is then condensed to butyric and caproic 
acids. 

The following information was obtained concerning the distribution of 
C4 in the fatty acids at the end of the fermentation. 


TABLE 3 
DISTRIBUTION OF C4 IN Butyric AcID 


Experiment 2 


PER CENT OF TOTAL 


CARBON ATOM .TS./MIN./MM. RECOVERED 
Carboxyl 1680 += 50 54 +2 
Alpha and gamma (average) 28 = 50 1+2 


Beta 1395 + 30 45 +2 


The final like the initial acetic acid contained C™ only in the carboxyl 
group. The actual results were 7400 + 200 cts./min./mM of carboxyl 
carbon and 0 + 30 cts./min./mM of methyl carbon. 

Data on the distribution of C' in butyric acid are presented in table 3. 
The butyric acid was obtained from a second experiment identical with 
experiment 1. The data show that Cis present in the carboxyl and beta 
positions but not in the alpha and gamma positions. The quantity of C' 
in the carboxyl and beta carbon atoms is almost equal, though there ap- 
pears to be a small but significantly greater amount in the carboxyl carbon. 
The smaller C'* content of the beta carbon may result from the fact that 
the beta (and gamma) carbon atoms are derived from an oxidation product 
of ethanol which is not in complete isotopic equilibrium with the preformed 
acetic acid. 

No method was available for determining the complete distribution of 
C'*in caproic acid. The presence and amount of C' in the carboxyl group 
could, however, be determined by decarboxylation of the barium salt. 
We found 4150 cts./min./mM in the barium carbonate obtained by the 
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decarboxylation of barium caproate containing 11,600 cts./min./mM. 
Therefore, 35.7 per cent or slightly more than one-third of the C'4 was pres- 
ent in the carboxyl group. The data are consistent with the view that the 
C* is equally distributed among three of the six carbon atoms (probably 
those in the carboxyl, beta and delta positions), but this is not certain. 

The results of experiment 1 (table 1), considered in conjunction with 
other information now available, lead us to postulate that the following 
reactions may be involved in the formation of butyric and caproic acids 
from acetic acid and ethanol. 


CH;CH.OH + K,HPO, — CH;COOPO;K, + 4H (8) 
CH;COOPO;K, + ADP = CH,COOH + ATP (9) 


In the presence of synthetic CH;C*OOH, the forward reaction (9) would 
cause a dilution of the active acetic acid and the reverse reaction would 
introduce C** into the carboxyl group of acetyl phosphate. Although there 
is at present no direct evidence to indicate the occurrence of acetyl phos- 
phate in this fermentation, the compound is almost certainly formed by 
Cl. butylicum,® an organism which also produces butyric acid (and butyl 
alcohol) from acetic acid. Moreover the chemical properties of acety] 
phosphate are consistent with its being an intermediate in fatty acid 
synthesis.”»* The labeled acetyl phosphate and the labeled acetic acid are 
assumed to condense to yield a product that is reduced to butyric acid. 
This reaction is not reversible. 


CH;C*OOPO;Ke + CH;C*OOH + 4H — 
CH;C*H,;CH,C*OOH + KzHPO, + H2O. (10) 


A similar type of condensation and reduction involving butyric acid and 
‘acetyl phosphate could lead to the formation of caproic acid labeled in the 
carboxyl, beta and delta positions. 


CH;C*OOPO;Ke + CH;C*HsCH;C*OOH + 4H > 
CH,C*H,CH,C*H:CH,C*OOH + KsHPO, + HO. (11) 


Reaction (11) represents only one of two possible ways in which caproic 
acid could be formed. Instead of there being a direct addition of the acetyl 
group to the terminal methyl group of butyric acid as is indicated above, 
the high energy phosphate group could be first transferred to butyric acid 
to give butyryl phosphate which could then condense with the methyl 
group of acetic acid. 

CH;C*H,CH2C*OOH + CH;C*OOPO;K; = : 

CH;C*H2CH2C*OOPO;K, + CH;C*OOH. (12) 


CH;C*H:CH,C*OOPO;Kz + CH;C*OOH + 
4H — CH;C*H,CH,C*H,CH2C*OOH + K2HPO, + H20. (13) 
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Evidence has been presented by Koepsell, et al.,° that a transphosphoryl- 
ation between acetylphosphate and butyric acid occurs under the influence 
of an enzyme preparation obtained from Cl. butylicum. With either 
mechanism the isotope distribution of the resulting caproic acid would be 
the same under the conditions of experiment 1. 

It is possible, however, to distinguish between these two types of con- 
densation by carrying out a fermentation of ordinary ethanol and synthetic 
butyric acid labeled in the carboxyl group. The position of the C' in the 
resulting caproic acid must depend on the way in which the condensation 
occurs. This is illustrated in equations (14) and (15) where the inactive 
acetyl phosphate and acetic acid are assumed to originate from the oxida- 


CH;COOPO;Kes a CH;CH.CH.C*OOH + 4H — 
CH;CH,CH:CH:CH2C*OOH + Ke,HPO, + H.O. (14) 


CH;(CH2)2C*OOPO;K2 + CH;COOH + 
4H — CH;(CH2)sC*H2,CH,COOH + K,HPO, + H20. (15) 


tion of ethanol. It can be seen that a condensation of acetyl phosphate 
and labeled butyric acid would yield carboxyl-labeled caproic acid, while a 
condensation of labeled butyryl phosphate and ordinary acetic acid would 
yield caproic acid labeled in the beta position. 


TABLE 4 
THE FERMENTATION OF CH;CH,OH anp SynTHETIC CH;CH,CH2C*OOH By 
Cl. kluyveri 
Experiment 3 
TOTAL 
COMPOUND MM/10 ML. cts./MIN./MM cTs./MIN, 
Before growth 
Acetic acid 0.100 0 0 
Butyric acid 0.384 17,200 += 500 6610 
Ethanol ca. Di 0 0 
After growth 
Acetic acid 0.101 0 + 100 0 
Butyric acid 0.05 18,000 = 10,000 955 
Caproic acid 0.427 12,600 = 500 5400 
Ethanol 0.306 250 = 50 75 


The results of an experiment designed to distinguish between these 
possible condensation mechanisms are given in table 4. Alcohol was 
added in excess to favor a high conversion of butyric acid to caproic acid. 

It can be seen in table 4 that most (81.7 per cent) of the C" originally 
present in butyric acid was converted into caproic acid. The molar 
activity of the caproic acid (12,600 + 500 cts./min./mM) is somewhat less 
than that of the original butyric acid (17,200 + 500 cts./min./mM) due 
to the fact that part of the caproic acid was completely built up from in- 
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active ethanol and acetic acid. This is evident from the fact that the 
quantity of caproic acid formed (0.427 mM) was greater than the butyric 
acid originally added (0.384 mM). When this dilution is taken into ac- 
count, the calculated molar activity of the final caproic acid (17,200 cts./ 
min./mM X 0.384 mM/0.427 = 13,400 cts./min./mM) agrees fairly well 
with its observed molar activity (12,600 cts./min./mM). 

The caproic acid formed in this experiment contained little or no C" in 
the carboxyl group. The barium carbonate obtained by decarboxylation 
from two samples of barium caproate showing activities of 845 + 20 and 
657 + 15 cts./min., were found to give only 6 + 3 and 5 + 3 cts./min., 
respectively. It will be recalled that synthetic carboxyl-labeled barium 
caproate yields 52 per cent of its C4 as barium carbonate. 

The absence of C' from the carboxyl group of caproic acid proves that a 
condensation of the type shown in equations (11) and (14) does not occur. 
The results are, however, consistent with the mechanism shown in equa- 
tions (13) and (15), involving a condensation of butyryl phosphate with 
the methyl group of acetic acid. The evidence in favor of this type of con- 
densation would be stronger if the C* had been proved to be present only . 
in the beta position. This has not yet been accomplished, but further 
work on this problem is in progress. 

There are a few other points concerning experiment 3 (table 4) that 
should be. mentioned. 

The complete absence of activity in acetic acid definitely proves that the 
conversion of acetic to butyric acid is not appreciably reversible at least so 
far as the carboxyl carbon atom is concerned. Evidence has already been 
presented above against any considerable exchange of carbon between 
acetic acid and the alpha, beta or gamma carbon atoms of butyric. 

A small but definite amount of activity was found in the neutral volatile 
compounds, reported as ethanol. It is possible that this activity is actually 
present in ethanol, though this would be difficult to reconcile with the fact 
that there is no activity in acetic acid which seems to be an intermediate 
between butyric acid and ethanol. Another possibility is that the activity 
attributed to ethanol is actually present in an associated compound such as 
butyl or hexyl alcohol. Since the molar activities of the C, and C. com- 
pounds are relatively very high, a small (1-2 per cent) contamination of 
the ethanol would account for the observed activity. With the small 
amount of alcohol and the low activity involved in the present experiment, 
the contamination theory could not be proved or disproved. However, 
it may be mentioned in this connection that the characteristic odor of C1. 
kluyveri cultures suggests the presence of a higher alcohol or ester. 

Summary.—1. When Cl. kluyveri is grown anaerobically in a medium 
containing ordinary ethanol and synthetic acetic acid labeled in the car- 
boxyl group with C™, labeled butyric and caproic acids are formed. 
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2. The butyric acid so formed has C** almost equally distributed be- 
tween the carboxyl and beta positions. The alpha and gamma positions 
are inactive. 

3. The caproic acid has one-third of its C'4in the carboxyl group; prob- 
ably the beta and delta positions are also labeled, though this has not been 
proved. 

4. No active carbon dioxide is formed from carboxyl-labeled acetic 
acid. This indicates that carbon dioxide is not an intermediate in these 
reactions. 

" §. The C1 content of the residual acetic acid is much lower than that of 
the initial acetic acid. This evidently results from the oxidation of ethanol 
to acetic acid or a related compound in isotopic equilibrium with it. 

6. When Cl. kluyvert is grown with ordinary ethanol and synthetic 
carboxyl-labeled butyric acid, C' is found in caproic acid but not in acetic 
acid. 

7. The active caproic acid so formed contains almost no activity in its 
carboxyl group. 

8. The above facts are consistent with the view that the formation of 
butyric acid involves a condensation between acetic acid or a reactive 
derivative thereof, such as acetylphosphate, formed by the anaerobic oxida- 
tion of ethanol, and a second molecule of acetic acid. The condensation 
product is then reduced to butyric acid. Caproic acid formation involves a 
condensation of the carboxyl group of butyric acid or some related C, 
compound, like butyrylphosphate, with the methyl group of acetic acid. 


* Present address: Mallinckrodt Institute of Radiology, Washington University, 
St. Louis, Mo. ; 
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THE DIFFERENTIAL EQUATION OF THE DISTRIBUTION OF 
GENE FREQUENCIES 


By SEWALL WRIGHT 
DEPARTMENT OF ZOOLOGY, THE UNIVERSITY OF CHICAGO 
Communicated October 15, 1945 


The first attempt to determine the mathematical form of the distribution 
of gene frequencies in populations was based on the setting up of differ- 
ential equations for certain special cases (Fisher, 1922,! 19307). A correc- 
tion and extension of these results came from expression of the conditions 
in an integral equation (Wright, 1929,* 1931‘). A general solution has 
since been obtained for fully stationary distributions by a third method 
(Wright, 1937,5 1938°). The case of uniform flux has been treated less 
generally (Wright, 1938,° 19427). Dr. A. Kolmogorov’ has recently been 
kind enough to send me a reprint of an important paper on this subject 
which was published in 1935 but which had not previously come to my 
attention. While the application is restricted to a particular stationary 
distribution, the method of approach points to a more systematic formula- 
tion than before. 

The situation discussed by Kolmogorov is that of a large population, 
consisting of many subgroups of size , each of which receives a certain 
number (k) of immigrants from the general population but otherwise 
breeds within itself. The average rate of change of the gene frequency of 
subgroups, in which is the frequency of a given gene, is represented by 
A = X(Ap) = (k/n)(p — p) where is the mean value of p in the whole 
population. The variance of p, due to accidents of sampling in one gener- 
ation, is represented by B = Z(Ap)? = pq/2n,q = 1 — p. It is stated, 
without demonstration, that the distribution u(p) of gene frequencies 
among subgroups after a stationary state has been reached, answers to the 
differential equation 


1 2 
2 dp? 





(Bu) — 5, (Aw) = 0. (1) 


The pertinent solution is given as 


u(p) = ptte—1gtka-1/B(4kp, 4kq). (2) 


The effect of selection in this situation is discussed briefly without, how- 
ever, modifying u(p) by introduction of the selection term, ap’q, into A. 

It is noted that the same formula (2) had previously been derived by the 
present author‘ by a different method. Equation (1) has, however, 
broader implications if valid for the general case A = 2(Ap), B = =(Ap)? 
and not merely for the particular case A = (k/n)(p — p), B = pq/2n. 
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The Immediate Factors of Evolutionary Change.—The immediate factors 
that tend to cause systematic changes (Aq) in gene frequency (gq) may be 
listed exhaustively* ® as (a) ‘mutation pressure, Ag = v(1 — gq) — uq, 
where v and w are the rates of mutation to and from the gene in question, 
(6) immigration pressure, Ag = m(q; — gq) (Kolmogorov’s A) where m is 
the proportion of replacement by immigrants and q; is the gene frequency . 
in these, and (c) selection pressure, which may take widely diverse forms 
but in the important case of constant relative selective values (W) for 
each multiple factor genotype in a random breeding population takes the 


OW , — Se : a asleht oa, 
form Aq = q(l1 — q) Oq/™: where 7 is 1 in haploids, 2 in diploids, the 


usual case, 1.5 for sex linked genes (if equal numbers of males and females), 
4 in tetraploids, etc.”.* In addition to these systematic pressures are (d) 
the random variations, 6g, due to accidents of sampling, the variance of 
which is Oa = g(1 — q)/rN in a population of effective size N (Kolmo- 
gorov’s B). The diploid case (r = 2) will be assumed in what follows. 

The Stationary Distribution of Gene Frequencies.—Systematic pressure 
toward the gene frequency, at which Ag = 0 and the cumulative effects of 
accidents of sampling determine a probability curve ¢(qg) describing the 
frequencies which would be exhibited in the long run by the value of q for 
a particular gene in a population subject to constant conditions. This 
distribution may also be interpreted as that exhibited at one time by the 
values of g in a group of populations that are all subject to the same condi- 
tions (as in the case of Kolmogorov’s u(p)). The deviations from the bi- 
nomial square formula for genotypic frequencies in the total population, 
depend on the variance of ‘y(g) under this interpretation.4:* '° In other 
cases ¢(q) may be used as the distribution at any time within either a class 
of non-allelic genes or an extensive series of multiple alleles,!! all subject to 
the same conditions. 

That equation (1) is, in fact, completely general for the stationary form 
of distribution may be shown by a slight modification of a method® ° that 
has been used for derivation of ¢(q). 

The conditions for stability of the distribution (including the terminal 
classes g = 0, g = 1) may be represented by two equations expressing the 
persistence of the mean and variance, respectively 


So'(q + 8g + Ag)e(g)dg = So'qv(q)d¢. (3) 
So'(a — q + 8¢ + Aq)*e(g)dq = So'(9 — 9)*¢(Q)da. (4) 


Noting that the mean value of 6g is zero, and that 6g is not correlated 
with g or Ag, these equations reduce to the following if the term in (Ag)? 
in (4) may be ignored. It may be noted in this connection that this term is 
negligible if Ag is of the same order as Ore or less, while if of higher order, 
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systematic pressure dominates the results so completely that the distribu- 
tion formula itself becomes unimportant. 


; So'Agqe(q)dq = 0. (3a) 
2 Si'(g — DAge(q)dg + So'o4,e(g)dq = 0. (4a) 

Putting Agy(q)dqg = dx(q) these conditions become 
x(1) — x(0) = 0. (30) 


2 So'x(q)dq — 2[x(1) + g(x(1) — x(0))] — So'e5,e(g)dq = 0. (46) 
Substituting (30) in (4b) the latter becomes 


So'[2x(q) — 2x(1) — o%e(g)]dq = 0. (4c) 


A solution is obtained by removing the integral sign since the resulting 
equation not only satisfies (4c) but also (3b) (noting that _ =Oifg=0 
or if g = 1, there being no sampling variance without alternatives in the 
sample). 

x(q) — x(1) = 1/205,0(q). (5) 


This can be solved for ¢(q) by differentiating the logarithm of the left- 
hand number and making the appropriate substitutions.°® 


(gq) = (C/o? eS aa/7oadae | 
where C is a constant such that fo'o(g)dq = 1. 


Since g increases by steps of 1/2N in a population of size N, the frequency 
of a given value of q is f(g) = ¢(q)/2N. From a study‘ of simple cases 
(N = 2 or 3) in which the frequencies in the stationary state can be deter- 
mined algebraically and from a more elaborate investigation by R. A. 
Fisher’ of the subterminal region in certain cases, it appears that the fre- 
quencies are given with considerable accuracy by the formula except for 
the terminal classes, g = 0,qg = 1. Consideration of the exchanges which 
occur between the terminal and neighboring classes leads‘ to the following 
approximate estimate for the terminal class, g = 0. That for g = 1 is 
analogous. 


f(0) = f(1/2N)/4N (mq, + »]. (7) 


The differential equation for the completely stationary case is given by 
differentiation of (5). It comes under equation (1). 


d 
qa (o309(q)) — Age(g) = 0. (8) 
q 


bole 


Since Agg(q) is the proportion of the distribution which tends to be 
carried past a specified value of g by the systematic pressure Aq, the other 
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term must represent the net proportion which tends to be carried in the 
opposite direction by accidents of sampling in each generation. 

The Case of Steady Flux.—There may be a practically stationary state of 
the proportions in all intermediate values of q¢ in spite of steadily increasing 
frequency of one terminal class at the expense of the other, provided that 
the proportion lost by the donor terminal class is negligible. This cannot 
be the case if either mutation rate or immigration rate is appreciable, but 
may hold in the presence of strong selection pressure since selection pressure 
is nil in populations in which g = 0 org = 1. 

The differential equation for the case of steady flux must differ from (8) 
by a constant term (D), the net proportion of the total (excluding the 
recipient class) that is carried past each value of g in each generation. 


d 
5 gg (otwe(@)) — Agel) + D = 0. @) 
q 


This is the general form given by one integration of (1) which is therefore 
the general differential equation for a steady state of the intermediate 
classes. It may be reduced to a linear equation of the first order by making 
the substitution, y = o3¢(q). 


d A 
< (2), + 2D =0. (10) 
dq oq 


The solution for ¢(g) is as follows: 


o(q) = [e 2S (49/05. 40/425) ic- 2D S e~*F 49/059 49q9) (11) 


The simplest special case is that in which Ag may be treated as zero 
(although there could be no flux if it were absolutely zero). 


C ra 2D 
Ging: trad 


The case under (12) that is most important genetically is that of irre- 
versible mutation at a rate so low that the donor class (¢ = 0, or g = 1) is 
not appreciably depleted. According to direction of mutation, 


f(Q) = 2v/q, or f(g) = 2v/(1 — q). (13) 


The ratio of the subterminal classes (1/2 in this case) gives the probabil- 
ity that a single neutral mutation may reach fixation instead of elimination. 

Returning to (12) the case in which D = 0 yields the corresponding sim- 
plest solution for a completely stationary state 





f@ = (12) 
q 


¢(q) = 1/[2(0.577 + log 2N)q(1 — q)] (terminal classes excluded). (14) 
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The case in which there are constant relative selection coefficients for all 
ow. — 
genotypes (Ag = g(1 — q) | rW) gives an apparently simple but in 
qd 


general rather refractory form (assuming a given set of frequencies of other 
genes) 


9(q) = [W°"/osd] [C — 2D W-?% dq]. (15) 


It will be convenient for later reference to cite the less general case 
Aq = q(1 — g)(s + éq), oc = g(1 — g)/2N which allows for any degree of 
dominance, provided s and ¢ are both small. 


fq = [eM /Q(1 — g][C— DS e~ENH2NHGg], (16) 


Non-stationary States —The general case, in which the proportion at 
each value of g is a function of time as well as of q itself, is given by the 
following, of which equation (1) is the case in which the left-hand member 
is zero. Time (7) is measured in generations. 


de(@, T) _ 1 2* 


2 ra) 
oT 2 dg? [osep(q, T)] — a [Age(q, T)]. (17) 


This can be reduced to an ordinary differential equation in the case in 
which the distribution has reached stability of form, with all classes (except 





the terminal ones) falling off at the same rate. Let K = — Me OE, 
9qT) of 
be the rate of decay per generation. 
be. % d 
5 7,2 (75e(Q)) — > (Age(g)) + Ke@ = 0. (18) 


2 dq’ dq 


It may easily be verified that for the case in which fixation is occurring 
under the uncomplicated effect of inbreeding (Ag = 0, K = 1/2) the only 
solution that does not involve negative frequencies is 


gq) = 1, orf(g, T) = Ce’, (19) 


In the case of irreversible mutation at an appreciable rate, Ag = v(1 — q) 
the rate of decay is easily shown to be K = v. Equation (18) is satisfied 
by the following value, originally derived by a different method. 


f(q@) = 2vqt?-?. (20) 


An analogous solution applies to the effect of swamping by immigration 
from a population in which the gene in question is fixed (Ag = m(1 — q), 
K = m) 


f(@) = 2mqt*™—", (21) 
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Comparison with Results by Other Methods.—The first attempt at deter- 
mining the distribution of gene frequencies was made by R. A. Fisher! 
who arrived at differential équations for certain special cases, in terms, 
however, of a different variable than gene frequency, ® = cos—! (1 — 2g), 
used in order to make the sampling variance constant. A discrepancy be- 
tween the rate of decay (K = 1/4N), derived by him for the case in which 
Aq = 0, and the value, 1/2N, given by a general method’ '* for determining 
the rate of fixation of genes under any system of mating, led the present 
author* ‘ to a different approach. The condition for a stationary state of 
the intermediate classes except for possible decay at rate K, was repre- 
sented by the following equation in which gq and x are recipient and donor 
classes, respectively, in the exchanges which occur from one generation to 
the next. 


00) _ 

2N 

(2N)! 
(2Nq)![2N(1 — g)]! 


It could easily be seen that if Ax = 0, the equation is satisfied by y(q) = 
g(x) = 1, K = 1/(2N + 1), the latter at least a close approximation to the 
rate of decay expected in this case. For the simplest stationary state, 
K = 0, dq = 0, the expression ¢(g) = Ag! + B(1 — q)~! is indicated 
(cf. 12). Approximate solutions could also readily be obtained for the 
linear pressures of mutation and migration. Selection presented more 
difficulty. 

On inspection of these results in manuscript, Fisher? was able to correct 
and extend his equations to obtain the following: 


(1 — K) 


Six + Ax)?%4(1 — x — Ax)?XO-VG(~)dx. (22) 





I. Case of uniform decay (Ag = 0) 


oy ligy . 0 | - 
Sans I os : 23 

oT = 4n E * 39 ate a 

y = Age 7" sind — (ef. (19)). (24) 


II. Stationary state, no selection (Ag = 0) 
> + ycot } = —4nB (25) 
do 7 


y = A cosec 3 + 4nB cot 3 (general, cf. (12)). (26) 
y = A cosec 3 (symmetrical case, cf. (14)). (27) 
y = 4nB(cosec 3 + cot #) (unidirectional mutation, cf. (13)). (28) 
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III. Stationary state, selection, no dominance, Ag = ag(1 — q) 


d 
A — (2an sin # — cot d)y = —4anA. (29) 


y = cosec #(2A + Be~20" 5 *) (general, cf. (16),¢ = 0). (80) 
(1 ss e 7 2an(1+cos 8) 





y = 4 cosec 3 (unidirectional mutation). (31) 


(1 = hee 

In cases I and II, the results agreed with those obtained from the integral 
equation (22), as may be seen by making the substitutions cos } = 1 — 2g, 
yd = 9(q)dq, d8/dqg = 1/~/q(1 — gq). In case III it was the author’s 
turn to make a correction in the selection term (published first® as e?%**), 
by taking cognizance of a series of small terms erroneously thought to be 
negligible but which actually doubled the exponent. With this correction, 
there was agreement.‘ 

The most general result® obtained for the completely stationary state by 
solution of (22) took into account all of the factors of change in the form 
Ag = v(1 — g) — ug— mq — gq) +911 — g)(s + 4), 6% = Q(1 — g)/2N. 


9(q) én Cpteet ee peat) yee ret etd" (32) 


This agrees with that obtained by substituting these values of Ag and of 
Ore in (6). 

The most general result®* obtained by this method for the case of steady 
flux was for Ag = g(1 — q)(s + ég). 


f(g = [eANet?N79(1 — g)|[C — 2Dge~ ON*FNY(2Nsq, 2Nig*)]. (33) 


where 
a a «@ — 
¥40)-1+ 5+ 7 t+7--=@-e )/2a 


b 7b? 2768 " 
¥(0, 6) = 1 re ee A 


Em = (Em—1 + Em—2)/2m(m + 1). 


No recurrence formula was recognized for the joint terms, ¥(2Nsq, 
2Ntg*) but the coefficients were calculated’ up to those pertaining to q’. 

The probability of fixation of a single mutation (C = 2v, D = ye NstMt } 
¥(2Ns, 2Nt) for irreversible mutations from class g = 0, or C= 0, D = 
—ve~ @NstN) 1Y(2Ns, 2Nt) for irreversible mutations from class qg = 1), 
could be calculated from the ratios of the subterminal classes, (Prob. = 
~/s/2N for a recessive mutation with selective advantage s, Prob. = 2s 


for a dominant mutation with selective advantage s, or for a semidominant 
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with selective advantage s in the heterozygote). The last agrees with 
Fisher’s conclusion.? Equation (33), with ¢ = 0, is indeed equivalent to 
(31). 

Comparison of (33) with (16) shows that if the former is correct, the 
following must hold: 


¥(2Nsq, 2Nig?) = (NF NI/g) Se“ GNAFINMdg. (34) 


This was tested by expanding the two exponentials in (34), integrating 
each term of the second one and combining. The coefficients were in all 
cases identical with those published’ for ¥(2Nsq, 2.Ntgq?). 

Equation (22) also gave the solution (20) for the case of uniform decay 
under an appreciable mutation rate. 

The integral equation (22) and the differential equation (18) are clearly 
equivalent to a close approximation. They are not exact mathematical 
equivalents, however, as may be seen from the fact that K must be put 
1/(2N + 1) in (22) if Ax = 0 to give the solution g(g) = 1, while it takes 
its true value 1/2N in (18) to give the same result. In the other cases (ex- 
cept (12)) second order terms have been omitted in the series, obtained as 
solutions of the integral equation, which do not appear in the solutions of 
the differential equation. Neither equation, of course, represents the natu- 
ral conditions exactly since integration is substituted for summation and 
differentials for minimal steps (1/2) in gene frequency. 


1 Fisher, R. A., Proc. Roy. Soc. Edinburgh, 42, 321-341 (1922). 
? Fisher, R. A., Zbid., 50, 205-220 (1930). 

3 Wright, S., Amer. Naturalist, 63, 556-561 (1929). 

4 Wright, S., Genetics, 16, 97-159 (1931). 

5 Wright, S., these PROCEEDINGS, 23, 307-320 (1937). 

6 Wright, S., Ibid., 24, 253-259 (1938). 

7 Wright, S., Bull. Amer. Math. Soc., 48, 223-246 (1942). 

8 Kolmogorov, A., C. R. del’Acad. des Sciences del’ U.R.S.S., 3 (7), 129-132 (1935). 
9 Wright, S., these PROCEEDINGS, 24, 372-377 (1938). 

10 Wahlund, S., Hereditas, 11, 65-106 (1928). 

11 Wright, S., Genetics, 24, 538-552 (1939). 

12 Wright, S., Ibid., 6, 111-178 (1921). 

13 Wright, S., Amer. Naturalist, 61, 330-338 (1922). 








390 PHYSIOLOGY: HALL AND BLAKESLEE Proc. N. A. S 


EFFECT OF SMOKING ON TASTE THRESHOLDS FOR 
PHENYL-THIO-CARBAMIDE (PTC) 


By Apa R. HALL AND ALBERT F. BLAKESLEE 


DEPARTMENT OF ZOOLOGY AND PHYSIOLOGY, WELLESLEY COLLEGE, AND DEPARTMENT 
oF BoTANy, SMITH COLLEGE 


Communicated November 2, 1945 


In the course of experiments on tasting ability the question arose as to 
whether a person’s threshold was the true one if he had been recently 
smoking. Salmon and Blakeslee (1935) tested a large group for PTC 
threshold and recorded the type of smoking for each individual (heavy, 
moderate, none). These records show that the position of the subject’s 
threshold on the PTC scale is not correlated with his smoking habits. 
That is, a heavy smoker is just as likely to have a low threshold as a high 
one and vice versa. The amount that an individual might vary from his 
resting threshold when he smoked again had not been tested. (Resting 
threshold is here used as the value obtained after nine or more hours of 
abstinence from tobacco.) 

Such a series of tests was therefore undertaken to show the effect, if 
any, that tobacco has on the taste apparatus, and to determine how long a 
time must elapse between smoking and a return to the original threshold. 

In a search of the literature for the effect of tobacco on the various senses 
it was found that the following reactions have been studied: visual acuity, 
skin pressure sense, eye accommodation, skin reaction to slight electric 
currents, taste preferences, and mental and physical efficiency. In addi- 
tion a number of workers have analyzed smoke and the tobacco in its 
various forms, snuff, pipe cuts, and cigarettes, for the ingredient causing 
the reactions. Bogen (1936) after careful analyses reports that no single 
explanation applies to the problem of irritation by tobacco smoke. The 
innocent bystander receives the sidestream rich in ammonia and other 
alkaline substances quite irritating to the eye and nasal passages. The 
smoker who gets the main stream through the cigarette or pipe stem has a 
more acid material thereby receiving his nicotine (alkaloid) as salts which 
ere less irritating. The smoker does get heat, pyridine, volatile acids, 
tarry and phenolic constituents, and the aldehydes, furfural and acrolein. 
These all lead to irritation of the membranes. All workers are fairly well 
agreed that nicotine in whatever form absorbed is the principal toxic agent 
wherever the nervous system is involved. 

Carefully controlled tests on mental and physical efficiency before and 
after smoking show a definite correlation between change in efficiency and 
smoking. Hull (1924) points out the following facts: 
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1. There is a large and uniform increase in the tremor of the hand 
lasting an hour and 23 minutes. 

2. There is marked and uniform stimulation of the heart which is still 
present an hour and 40 minutes after smoking. 

3. There is a minute increase in the speed of reading reaction time. 

4. There is a gain in the rate of complex mental addition lasting an 
hour and 15 minutes but with no measurable effect on accuracy. 

5. There is a high probability of loss in auditory memory and efficiency 
in rote learning, with recovery in an hour. 


In a series of ten mental tests Bush (1914) reports a 10.5% decrease in 
efficiency after smoking. 

Accommbdation time, both near-to-far and far-to-near sight, has been 
studied over a period of years at Wellesley College. Homewood and 
Howe (1937) have shown that accommodation time is definitely decreased 
during the first 20 minutes after smoking, but increased from 40 to 60 
minutes after. This agrees with Schrumpf-Pierron’s (1927) statement that 
the effect of tobacco on the nervous system is first a stimulation and then 
a depression. They also found that the occasional smoker was stimulated 
more than the habitual one and depressed sooner and to a greater extent. 

Wenusch and Schiller (1936) working on skin pressure found that sensi- 
tivity to both hair pressure and pendulum stroke was changed during 
smoking. They do not record any stimulating effect of smoking, only a 
depressant one. 

From a survey of the literature on the changes in mental and physical 
states after smoking, and from a series of electrical tests on the finger before 
and after smoking, Mendenhall (1930) concludes that tobacco has the same 
effect as rest. If a person is tired and depressed, smoking will stimulate 
and bring him back to normal, if overexcited it will quiet him by its de- 
pressing action. If he is in a well-rested state it has no marked effect. 

Sinnot and Rauth (1937) found the thresholds for sugar and salt high 
in smokers, but during a period of several days during which six smokers 
had abstained, their thresholds fell to the level of non-smokers. Laird 
(1939) tested the effect of the smoking habit on the preference for sweet or 
tart pineapple juice. He found no difference in the taste preferences of 
smokers and non-smokers among men, at any age, nor among women up 
to forty years. Among the women of fifty to sixty-eight years the non- 
smokers were like the other groups, but the smokers were predominantly 
in favor of the tart juice rather than the sweet. He did not test each 
individual before and after smoking, nor does he deal with threshold 
ability. 

In testing the present group of subjects for the effect of smoking on 
tasting ability each person was his own control. Smokers were tested 
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only after nine or more hours of abstinence from tobacco. Of the 60 
subjects, 32 were from the staff at Cold Spring Harbor and 28 were ad- 
vanced students and staff at Wellesley College; 32 were habitual smokers 
and 28 non-smokers; 24 were men and 36 women. Concentrations of 
PTC were used from 1:5120 M up to 1:312.5 using a factor of two, such 
that each solution used in a test was twice as concentrated as that pre- 
viously administered. (M in these numbers represents 1000.) In each 
test approximately 0.6 cc. was given by means of the straw method. 


Subj. Gone. BeforelAfter | Min. after pmoking: Type of 
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FIG.|, ABILITY TO TASTE RTC. AFTER SMOKING. 





Two types of experiments were performed. The procedure for the first 
series of tests was as follows: 


1. The PTC threshold was determined. 

2. The subject smoked two cigarettes of a standard brand in 10-15 
minutes. 

3. The PTC threshold was determined immediately after smoking and 








ch 
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at 15-minute intervals for one to two hours. (In several cases the tests 
ran still longer.) In the second series 10 individuals of the first series, 
both smokers and non-smokers, were tested as above except that the smoke 
was drawn through a dry flask for cooling and into the nose through a steri- 
lized nose-piece keeping the mouth closed so that the smoke did not touch 
the taste buds. 


A summary of the results of the first series may be found in table 1. 
Figure 1 shows typical PTC curves for each of the behavior groups. 


TABLE 1 


SUMMARY OF THE REACTIONS TOWARD PHENYL-THIO-CARBAMIDE AFTER SMOKING 
(SMOKE TAKEN IN THROUGH THE Movuts#) 





INDIV. % OF NO. PTC RETURN TO INITIAL 
CHANGE IN TASTING IN TOTAL -—GRADES— LEVEL————~. -—SMOKER-— -——-SEX— 
ABILITY GROUP TESTED INC. DEC. NO. TIME YES NO fous Q 
A. None + 6.6 §... «.. -@ aigehauge 2.3 : ae 
B. Decreased only 32 63.38  ... 1-5 22 30’-60’ io it * Ia" 24 
10 Over 60’ 
C. Increased only 10 16.7 1-4 6 15-60’ = 2 ae 
4 Over 60’ 
D. Decreasedthen 12 20.0 1-2 1-2 7 30’-60’ SF Gio7 
increased 5 Over 50’ 
E. Increased then 2 Sca. 24 1-2 2 Over 60’ > eee ee 
decreased 
Totals 60 100 35 15’-60' 32 28 24 36 
21 Over 60’ 


There were three types of reaction: 


1. No change in tasting ability, group A (6.6%). 

2. Decrease in tasting ability, groups B and D (73.3%): group B re- 
turned to initial level (53.3%); group D returned to initial level and then 
showed an increase (20%). 

3. Increase in testing ability, groups C and E (20%): group C returned 
to initial level (16.7%); group E returned to initial level and then showed 
a decrease (3.3%). 


These figures check very well with those of Mendenhall (1930) working 
on sensitivity to electric shock. He found that when tested before and 
after smoking (750 observations) 72.2% of the cases showed a depression 
while 28.3% showed a stimulation. In our work for taste 73.3% were 
depressed, 20% were stimulated and 6.6% showed no change. 

Considering the time at which depression and stimulation occur the 
following summary may be made: 
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Time at which tasting ability is first lowered (groups B, D): 


38 during smoking (10-15 minutes) 
4 by 15 minutes after smoking 
2 by 30 minutes after smoking 


Time at which tasting ability is first increased (groups C, E): 


8 during smoking 
3 by 15 minutes after smoking 
1 by 30 minutes after smoking 


Time of greatest increase after an initial decrease (groups B, D): 


5 by 15 minutes after smoking 
10 by 30 minutes after smoking 
16 by 45 minutes after smoking 
10 by 60 minutes after smoking 
3 not increased again by 60 minutes (one was still low at 105 min.) 


Thus it may be seen that the initial depression occurs during smoking or 
by 15 minutes after. For those who show no initial decrease, stimulation 
occurs for most individuals during smoking or in 15-30 minutes after. If 
there is an initial decrease the return to the original threshold or to a period 
of stimulation is for most subjects at 30-45 minutes after smoking. 

The taste apparatus according to Ranson (1939) is composed of taste 
buds on the tongue from which two sets of special afferent visceral fibers 
pass to the tractus solitarius via the chorda tympani (of the seventh cranial 
nerve) and the glosso-pharyngeal nerve. These taste fibers connect with 
the anterior part of the nucleus of the tractus solitarius. Further con- 
nections are rather vague but certainly reflexes to the gustatory center and 
the motor nuclei for mastication and swallowing exist. A cerebral center 
is not definitely fixed but a spot near the anterior end of the temporal lobe 
is thought to have such function. There are at least three synapses then 
between the surface of the tongue and consciousness. 

According to the work done on other sensory-motor arcs the expected 
curve for taste should show an initial level, a more sensitive period, and 
then a depression period. But we have here an added factor to complicate 
results, namely, the exposure of the sensory endings, the taste buds, to the 
various drugs resulting from the combustion. It may be that the de- 
pressant effect on the taste buds antedates or coincides with the stimulative 
effect on the nerve cells and thus we get several forms of curves for the first 
half hour depending on the comparative strength of the two factors. 
Some individuals keep to the same level at first with later depression, 
others have a sharp early depression followed by stimulation, while others 
are stimulated immediately and to a greater degree than the depression 
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of the taste buds. We felt that it might be enlightening so to conduct the 
smoking period that the taste buds would not be touched by the smoke. 
Of the 10 individuals so tested 3 showed no change in tasting ability for 30 
minutes to an hour followed then by a depression. Six had an initial 
stimulation followed at 30 to 45 minutes by a depression. These persons 
had all shown initial depression when smoking by mouth. One subject 
(of group C) who showed only stimulation when smoking by mouth had 
the same type of reaction (stimulation only) when smoking by “nose.” 
Figure 2 shows the mouth and “nose” curves for two of these subjects. 


Subj- Core. a Min. after smoking: 

No. smokecig.| 15 | 30 |45 {60 | 75 | 105 | 120 
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FIG.2. EFFECT OF SMOKING BY NOSE. 


It would seem then that the sharp depression noted for the larger number 
of subjects in series 1 may be due to the direct action of some product of 
combustion on the taste buds rather than to the effect of nicotine on the 
nerve cells. 

We may therefore make the following conclusions: 


1. There is a definite effect on the individual’s ability to taste PTC 
after smoking, 73.3% of the subjects requiring stronger solutions in order 
to taste after smoking and 20% of them tasting weaker solutions. 

2. The time which must elapse after smoking before the individual is 
a proper subject for tasting experiments varies with the individual. In 
the present series of tests only 58% had returned to the resting threshold 
within an hour, while some took several hours for recovery. 

3. The initial effect of smoking in the larger number of cases is a direct 
dulling of the taste buds by some product of the combustion. 
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4. The true effect of nicotine on the nerve apparatus for taste appears 
to be the same as for other nerves tested—an initial stimulation with later 
depression. 
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RESPIRATORY ENZYMES IN PARAMECIUM: 1. CYTOCHROME 
OXIDASE 


By Epcar J. BoELL 
OSBORN ZOOLOGICAL LABORATORY, YALE UNIVERSITY 


Communicated October 15, 1945 


It has been generally stated, as the result of the work of a number of 
investigators,': *» * that cyanide is without effect on the respiratory activity 
of Paramecium. The report of Kalmus‘ that respiration of Paramecium 
is depressed by cyanide has been largely disregarded because of serious 
defects in his experimental technique (cf. Howland and Bernstein‘). 
Cyanide was likewise shown to be without effect on the respiration of a _ 
number of other ciliates,®: 7 * and the conclusion was reached, as sum- 
marized by Lwoff,* that “Cette insensibilité 4 HCN et a CO n’est pas 
générale chez les Protozoaires, mais est jusqu’ici particuliére aux In- 
fusoires.”” The repeated failure to obtain depression of respiration in 
ciliates with cyanide and similar inhibitors led naturally to the belief that 
oxidations in these forms were mediated by a mechanism which was 
different from that of most aerobic cells.’ 
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More recently, a reinvestigation of the effect of cyanide on the respira- 
tion of certain ciliates has been made. In these studies a number of 
technical improvements, not alone in the use of cyanide but also in rearing 
the animals and preparing them for respiratory measurements, were 
employed. Thus Baker and Baumberger™ and Hall'! showed that res- 
piration in Tetrahymena gelew and Colpidium campylum could be inhibited 
by cyanide. Kitching,!* in a study of contractile vacuole activity, ob- 
tained indirect evidence for cyanide susceptibility in the peritrich, Zo- 
othamnion. Finally, the respiration of Paramecium calkinsi was shown to 
be inhibited to approximately 50 per cent of the normal by cyanide (Boell!* 
and recently confirmed by Pace") and by azide of the same concentrations 
as those usually employed for respiratory studies. 

Since these inhibitors have been widely used to indicate the operation 
of the so-called cytochrome-cytochrome oxidase system in normal re- 
spiration, the conclusion was drawn" that in so far as susceptibility to 
cyanide and azide can be used as tests, the results suggest that the re- 
spiratory mechanism of Paramecium resembles that of most aerobic plant 
and animal cells. The purpose of this paper is to provide further support 
for the above conclusion by presenting direct evidence for the occurrence 
in Paramecium of an oxidative enzyme which presumably is cytochrome 
oxidase. Apparently this represents the first time that an oxidative 
enzyme has been obtained from a eiliate and studied by quantitative 
techniques. 

Technique.—Respiratory measurements and cytochrome oxidase de- 
terminations were carried out in the Cartesian diver apparatus on known 
numbers of Paramecium calkinsi which had been reared and prepared for 
study by methods previously described by Boell and Woodruff. In the 
cytochrome oxidase determinations, highly concentrated suspensions of 
well-washed Paramecium were used. Duplicate 10-cu. mm. samples were 
withdrawn in order to obtain a count of the number of animals present in 
a given volume of medium, and the remainder of the suspension (approxi- 
mately 0.1 cc.) was ground into a cell-free brei by means of a motor-driven 
homogenizer. Appropriate amounts of the enzyme preparation were then 
added to measured quantities of cytochrome-c, ascorbic acid (Merck’s 
Cebione) which had been neutralized with NaOH, and phosphate buffer 
at pH 7.2. A series of divers was then prepared, each containing a 1.5-cu. 
mm. sample of the mixture or the same volume of a preparation containing 
enzyme which had been heated in a boiling water bath for 5 to 10 minutes 
prior to mixing with cytochrome and ascorbic acid. In some cases divers 
were prepared so that comparisons of enzyme activity in the presence and 
absence of added cytochrome could be made. All experiments were 
conducted at 25°C. ‘ 

Results —1. Inhibition of Respiration by Cyanide*and Azide.—The 
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inhibition of respiration of intact cells by cyanide, azide, carbon monoxide, 
etc., has generally been interpreted as due to interference with the func- 


tioning of cytochrome and cytochrome oxidase. 


Thus, susceptibility to 


cyanide has been used as an indicator of the operation in vivo of this system 


of respiratory catalysts. 


The results of a typical experiment with these 
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inhibitors are shown graphically in figure 1. Curves 2 and 4 of this figure 
show, respectively, the effect of sodium azide and potassium cyanide, each 
in a final concentration of 10-? molar, on the oxygen consumption of 
Paramecium. With both these substances respiration is depressed to 
approximately 50 per cent of the normal, and data are available which show 
that this inhibition is completely reversible. At acid pH values, however, 
Paramecia are irreversibly affected by strong concentrations of azide. 
Such animals apparently lose their ability to maintain water balance, for 
the contractile vacuoles stop beating and become enormously swollen. 

In some experiments in which azide was used, there was a period of 
accelerated oxygen consumption following the inhibition of respiration. 
In Curve 2, representing only the first 75 minutes of a longer experiment, 
such a rise in respiratory rate can be 20 
seen. The rate of respiration stimu- 
lated by azide may, under certain 
experimental conditions, consider- 
ably exceed the normal, and it is 
of interest to note that the accelera- B. 
tion of oxygen consumption pro- a’ => rn > 
duced by azide is susceptible to oy 30 45 
epuaiibe. MINUTES 


FIGURE 2 
The control respiratory rates, as Oxidation of p-phenylenediamine by 
shown by Curves 1 and 2in figure 1, oxidase from Paramecium calkinsi. Ex- 
although identical, were obtained perimental conditions: 25°C.; pH 7.2; 
from two different samples of Para- mal concentrations of cytochrome-c = 
‘ a 1 X 107‘ molar, of p-phenylenediamine = 
mecium and at different pH values, 0.033 molar. Oxidase preparation equiva- 
and, incidentally, it may be men- ent to approximately 75 Paramecia. 
tioned that the average respiratory Curve A = oxidase + cytochrome + 
rate for these animals is essentially -phenylenediamine; Curve B = same 
identical with that reported for P. but with boiled oxidase. Ordinate denotes 
calkinsi of the same race and *8*? uptake in mal. 
mating type by Boell and Woodruff." 

2. Cytochrome Oxidase Activity.—Preliminary tests of the cytochrome 
oxidase activity of Paramecium were performed with 0.033 molar p-phenyl- 
ene diamine as substrate (Fig. 2). However, this compound, although it 
was readily oxidized, appeared to be toxic to the oxidase as evidenced by 
a rapid decline in the rate of O2 uptake (cf. Schneider and Potter'®). Its 
use was therefore abandoned. Ascorbic acid in a final concentration of 
0.01 or 0.02 molar was used instead, since a number of workers had shown 
that it can reduce cytochrome-c and could, accordingly, be used as a sub- 
strate in tests for cytochrome oxidase.'*— 

The results of a typical experiment are presented graphically in figure 3 
and show unmistakably the presence of an enzyme in Paramecium capable 
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cf oxidizing ascorbic acid. Moreover, since the catalytic oxidation of 
ascorbic acid is much greater in the presence of cytochrome than in its 
absence, it seems safe to conclude that the enzyme in this case is cyto- 
chrome oxidase. In this connection it.may be recalled that Sato and 
Tamiya™ reported finding the absorption bands of reduced cytochrome in 
Paramecium. 
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MINUTES 
FIGURE 3 
Oxidation of ascorbic acid by oxidase from Paramecium 
calkinsi. Experimental conditions: the same as stated for 
figure 2 except that ascorbic acid in a final concentration of 0.01 
molar was used as substrate in place of p-phenylenediamine. 
Curve A = oxidase + cytochrome + ascorbic acid; Curve 
B = same but with boiled oxidase; Curve C = oxidase + 
ascorbic acid; Curve D = oxidase + cytochrome; Curve E = 
oxidase. Ordinate denotes oxygen uptake in mul. 


Ascorbic acid in solution shows a fairly rapid oxygen uptake even in the 
absence of enzyme which has been shown to be due to the presence in the 
reaction mixture of metallic copper impurities.» °° It is necessary, there- 
fore, to run a blank containing boiled enzyme, cytochrome and substrate 
in order to determine the non-enzymic oxidation of ascorbic acid. How- 
ever, since Stotz, et al., have shown that the catalytic effect of metallic 
copper is eliminated by the presence in the reaction mixture of undenatured 
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protein, as well as other organic compounds,”! e.g., glycine,”* it seems fairly 
certain that this method introduces an over-correction. Consequently 
the cytochrome oxidase values obtained are minimal. 

Schneider and Potter’ have described a method of correction which 
involves making determinations of the oxidase activity at three different 
concentrations and then extrapolating these to zero in order to determine 
the oxygen uptake due to non-enzymic catalysis. Such a procedure has 
been used in the preparation of table 1 which shows, in addition, a direct 
relationship between enzyme activity and enzyme concentration. 


TABLE 1 
NUMBER OF ANIMALS 
IN 1.5 cu. MM. OF O2 UPTAKE CORRECTED CYTOCHROME RATIO: 
REACTION MIXTURE myl/HOUR OXIDASE ACTIVITY ACTIVITY/NUMBER 
0 35* 0 ee 
102 57 22 21 
204 76 41 20 
306 96 61 20 


* Derived by extrapolating the O. uptake at other values to zero enzyme concen- 
tration.'* 


Summary.—By a micromanometric method, the respiration of Para- 
mecium calkinsi has been found to be sensitive to azide and cyanide. In 
addition, an enzyme capable of oxidizing ascorbic acid and p-phenylene- 
diamine in the presence of cytochrome-c and which may, therefore, be 
designated as cytochrome oxidase has been demonstrated. 
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OFFSPRING FROM UNBORN MOTHERS* 
By W. L. RUSSELL AHD PaTrRicIA M. DouGLass 
Roscoe B. JACKSON MEMORIAL LABORATORY, BAR HARBOR, MAINE 
Communicated October 29, 1945 


The preliminary experiments reported here show that mouse ovaries 
transplanted from embryos to adult females can develop to maturity and 
produce offspring. 

Eight such operations were performed, three of which have proved suc- 
cessful. The donor ovaries were obtained from embryos aged 13!/, days 
(in 3 of the operations) and 14!/2 days (5 operations) post coitus, matings 
being timed by the vaginal plug method. They were transplanted to the 
ovarian capsules of young adult animals whose ovaries had been removed 
immediately before. A male was placed with each host a few days after 
the operation. Other transplantation experiments have shown that 
ovarian regeneration frequently occurs in host animals.! In order to dis- 
tinguish offspring of the transplanted ovaries from offspring of regenerated 
ovaries a mating plan was used which has been described by Russell and 
Hurst.!. This plan uses an inbred strain, stock 129, which carries both 
chinchilla, c”, and albino, c’, genes. The embryo donors were obtained 
from matings of stock 129 albinos and were, therefore, of the genotype 
c’c’. The hosts were hybrids, although for the purpose of this experiment, 
transplantations couid have been made within the 129 stock. The hy- 
brids were the offspring of matings of 129 cc X C57 black strain and, 
therefore, had the genetic constitution Cc". These Cc" hosts with c*c* 
ovary implants were mated with 129 cc males. Offspring of regenerated 
host ovaries would, therefore, be either Cc“ or c“"c**, while offspring from 
the transplanted ovaries would be cc’. 

The animals used have not yet reached the end of their breeding period, 
but the three successful operations have so far yielded, respectively, 3, 4 
and 10 offspring from the transplanted ovaries. In these operations the 
donor ovaries were obtained from embryos aged 14!/, days pe and were 
transplanted to hosts 40 days old. Regeneration of host ovaries occurred 
to some extent in all three cases as well as in four of the unsuccessful 
operations. The remaining unsuccessful operation produced no offspring. 
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The demonstration that offspring can be obtained from transplanted 
embryonic ovaries, that is, from unborn mothers, introduces a method that 
should prove of value in several different fields of experimental zodlogy. 

_Some of the possible uses are discussed below. 

1. Several dominant gene mutations apparently produce, in the 
homozygous mutant type, a lethal condition causing death of the animal 
while still within the uterus of the mother. The evidence rests on the pro- 
portion and nature of the abnormal embryos produced by matings of the 
heterozygotes. Transplantation of ovaries of the lethal type should, first 
of all, prove whether the ovary can survive in a normal host and, second, if 
the ovary does function, it would make possible a direct breeding test of 
the genetics of the lethal type itself. One experiment of this kind is being 
made by the authors. 

2. In the study of chromosome translocations, experiments similar to 
the above could be undertaken to determine the chromosome make-up of 
embryos that die presumably from chromosome unbalance. 

3. Ovarian transplantation could be used to identify the genetic type 
of embryos that have to be killed for study. An example will perhaps 
make this clearer. Mice homozygous for the gene W die shortly after birth 
from severe anemia and can, therefore, be obtained only from matings of 
the heterozygotes. Griineberg* points out that since the anemia of WW 
embryos is already well marked on the\16th day it must have started earlier. 
Its actual onset has not been determined because there was no way of 
identifying the homozygous segregants before that age. 

4. Many transplantation experiments that were possible with older 
ovaries could employ embryonic ovaries to obtain additional information. 
For example, in studies of differences in maternal environment the foster 
mother’s influence could be brought into play at an earlier stage in the 
development of the ova. 

5. The relation of hormones and other factors to the development and 
functioning of the ovary could be investigated by transplanting embryonic 
ovaries to hosts of various ages or to hqasts subjected to various experimental 
treatments. 

6. It is possible that transplantation would prove successful with still 
earlier embryonic stages. In that case it might be a useful technique for 
investigations on the origin of the germ cells. 

7. There is no reason to suppose that the method could not be repeated 
for any number of successive generations, thus leading to an indefinite 
number of unborn direct female ancestors. This possibility should be of 
interest to investigators in several fields of research. In studies of the 
mammary tumor agent, for example, mice could be obtained whose female 
ancestors, for any number of generations, had not been nursed. 

Summary.—Offspring were obtained from transplanted embryonic 
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mouse ovaries. This introduces a method that could be used for the gen- 
etic identification of some lethals and certain types of embryos that have to 
be killed for examination. It should prove of value in studies on differences 
in maternal environments, in experiments on the relation of hormonal and 
other factors to the development and functioning of the ovary and possibly 
for investigations on the origin of the germ cells. An extension of the 
method over successive generations, to produce animals descended from 
any number of unborn direct female ancestors, would provide still further 
possibilities for research. 


* Aided by a grant from the Rockefeller Foundation and by grants to the Roscoe B. 
Jackson Memorial Laboratory from the National Cancer Institute, the Jane Coffin 
Childs Memorial Fund for Medical Research, the Anna Fuller Fund and the Inter- 
national Cancer Research Foundation. 
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INDUCTION OF CONJUGATION IN PARAMECIUM, BURSARIA 
AMONG ANIMALS OF ONE MATING TYPE BY FLUID FROM 
ANOTHER MATING TYPE* 


By TzeE-TuaNn CHEN 
DEPARTMENT OF ZOGLOGY, UNIVERSITY OF CALIFORNIA, Los ANGELES 


Communicated October 25, 1945 


In these experiments, it was found that animal-free fluid from a Russian 
clone (Ru22) induces clotting and conjugation among animals of another 
mating type, even though the latter belongs to a different variety.! 

Most of the present work was done on the effects of fluid from cultures 
of a Russian clone (Ru22) on animals of mating types belonging to four 
other varieties. This Russian clone (together with a number of other 
Russian clones) was sent to the present writer from Moscow by Professor 
G. F. Gause. Jennings and Opitz? found that this clone did not mate with 
any of the four varieties known at that time and they were of the opinion 
that this clone possibly belongs to a fifth variety. In the present paper it 
is called the fifth variety. (Apparently only one mating type of this 
variety has been found.) 

The animals of this Russian clone (as well as those of other clones re- 
ferred to in this paper) were cultured in essentially the manner described 
by Jennings.* 

In testing the effect of fluid from the Russian clone on animals of another 
mating type, the fluid is first taken out of the culture (Ru22) with a micro- 
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pipette; it is then placed in depression slides and carefully examined under 
a dissecting microscope. All the animals, if any are accidentally included 
in the fluid, are removed. Then a number of animals of the clone to be 
tested, together with a very small quantity of fluid from this clone, are 
added to the Ru22 fluid. The depression slides are then left in a moist 
chamber and examined from time to time. As a control, fluid from the 
clone to be tested is first taken out of the culture, placed in a depression 
slide, and then animals from the same culture (or another culture of the 
same clone) are added. 


Se 


ja 
FIGURE la 


Clotting and pair-formation among animals of an English clone (En1) 
induced by fluid from a Russian clone (Ru22). 


The fluid from cultures of the Russian clone (Ru22) was tested on ani- 
mals of a number of mating types belonging to four other varieties (II, III, 
IV, VI) and especially on animals of an English clone (En1) belonging to 
the sixth variety. 

Effects of the Ru22 Fluid on Animals of the English Clone (En1).—This 
English clone (Enly was collected by Professor E. G. Pringsheim in Cam- 
bridge, England, and was sent by him to the present writer, who subse- 
quently made numerous tests on it. The results so far obtained show that 
it will not mate with animals of any of the known varieties but will con- 
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jugate readily with animals of a clone collected from Prague, Czecho- 
slovakia (also sent by Professor Pringsheim). It seems clear that this 
English clone together with the Czechoslovakian clone constitute a new 
variety, which is designated as Variety VI. 

Clotting: The animal-free fluid from Ru22 cultures apparently renders 
the surface of the animals of the English clone sticky so that the animals 
adhere to one another, forming clots and pairs (Fig. la). No such clotting 
or pair-formation was ever observed in the controls. The clotting or group- 
formation of animals of the English clone (one mating type)‘ induced by 
the fluid of the Russian clone (another mating type) differs somewhat from 
the clotting which occurs when animals of diverse mating types are mixed. 


\ 
gZ 


Pairs formed among animals of an English clone (Enl) induced by the 
fluid from a Russian clone (Ru22). 


FIGURE 16 


The difference is twofold: (1) The clotting among animals of the same 
mating type occurs much later than that which usually occurs when ani- 
mals of diverse mating types are mixed. In the former case, agglutination 
occurs 4-6 (usually more) hours after the animals of one mating type are 
introduced into the fluid of another mating type; in the latter case, 
agglutination occurs immediately or almost immediately. (2) The clots 
formed by animals of the same mating type are much smaller than those 
formed as a result of mixing animals of diverse mating types. In induced 
group-formation, each clot usually consists of 2-7 individuals. Clots 
formed as a result of mixing animals of diverse mating types are usually 
larger; each of the larger clots often consists of many individuals. 
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Conjugation: Pair-formation among animals of the same mating type® 
differs somewhat from pair-formation when those of diverse mating types 
are mixed. In the first place, in induced conjugation, the pairs are formed 
much later (4-6 or more hours after the animals are introduced into the 
fluid of another mating type). When animals of diverse mating types are 
mixed, usually many pairs are formed in one or two hours. Secondly, in 
induced conjugation many of the pairs are atypical. Although some of the 
pairs are hardly distinguishable from those formed as a result of mixing 
animals of diverse mating types, other pairs are atypical in that the two 
conjugants do not have the relative positions typical of normal conjugating 
pairs (Fig. 1). Thirdly, in many pairs, the conjugants are not held as 
firmly together as those formed in ordinary conjugation. Asa result, ina 
number of cases the conjugants become separated when fixed. This is true 
even after the animals have been conjugating for more than 20 hours. 


aie 
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FIGURE le 


Single animals of an English clone (En1) undergoing 
perhaps autogamy or endomixis or some other type of 
nuclear processes induced by the fluid from a Russian 
clone (Ru22). 


Fixed and stained preparations of induced pairs show that nuclear 
changes occur in the conjugants and anlagen are formed in the ex-conju- 
gants. The nuclear changes in induced conjugation are being studied and 
will be reported on separately. 

Induction of Nuclear Changes in Single Animals: The Ru22 fluid 
induces not only clotting and conjugation, as already described, but also 
nuclear changes in some individuals not associated with any others. These 
single animals develop a paroral cone’(Fig. lc). This structure has been 
reported in Paramecium aurelia undergoing autogamy (endomixis).6 Pre- 
liminary cytological studies show that nuclear changes take place in these 
single animals. It remains to be determined whether they undergo autog- 








408 ZOOLOGY: T. T. CHEN Proc. N. A. S. 


amy, or endomixis, or some other type of nuclear processes. These single 
animals may be found together with the conjugating pairs; at other times, 
however, they may be found alone. 

In addition to the induction of clotting and conjugation among animals 
of the same mating type and of nuclear changes in some single animals, 
there are other effects of the Ru22 fluid on the members of the English 
clone. These additional effects are recognizable even before clotting and 
conjugation occur. The animals become (1) sluggish in movement, (2) 
darker in color, and (3) distorted in form to a greater or less extent. No 
such changes were found among the animals in the controls. 


Wea 
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FIGURE 2¢ 


Clotting and pair-formation among animals of an American 
clone, MM17 (Variety II), induced by fluid from a Russian 
clone (Ru22). 


It should be noted that only a percentage (usually small) of the animals 
are visibly affected by the Ru22 fluid, and hence only a small number of 
them undergo clotting and conjugation; while a few solitary animals 
undergo nuclear changes. It is difficult to understand why only a percent- 
age of animals are obviously affected. It seems probable that there are 
differences in physiological condition even among the animals of the same 
culture. It thus seems that those individuals which are affected differ 
from others added to the fluid at the same time. 
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Effect of Washing Animals.—In order to determine whether or not induc- 
tion of clotting and pair-formation among members of the English clone 
Enl by the Ru22 fluid might be due to some impurities in the cultures of 
the Russian clone, on November 23, 1944, a number of the animals of this 
Russian clone were carefully washed and then recultured. Some weeks 
later the fluid was tested again on the animals belonging to the English 
clone Enl. Clotting and pair-formation were again observed, as described 
above. 

Effects of the Ru22 Fluid on Animals Belonging to Some Additional 
Varieties—The fluid of the Russian clone Ru22 is apparently effective in 
inducing conjugation throughout the species since it induces conjugation in 
all the four varieties so far tested. The one other variety (I) is being 
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FIGURE 2b 


2b 


Pairs formed among animals of an American clone, MM17 
(Variety II), induced by the fluid from a Russian clone (Ru22). 


In addition to inducing conjugation among the animals of the English 
clone (Enl) belonging to the sixth variety, this Ru22 fluid also induces 
clotting and conjugation in each of the following clones: MM17 (mating 
type F) and Gr13 (type K) of the second variety, clones Grl (type N) and 
Pi3 (type O) of the third variety, and clone Ru3 (type R) of the fourth 
variety.’ Clotting and pair-formation in these varieties are essentially the 
same as those which occur among animals of the English clone as described 
above (cf. Figs. la, 16 and 2a, 2b). No clotting or conjugation was ever 
observed in the controls. - Some other clones of these three varieties were 
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also tested but as yet no clotting or conjugation has been observed. These 
negative results might be due to differential effects of the Ru22 fluid on 
different clones of the same variety. 

Discussion.—The phenomena in P. bursaria described in this paper are 
possibly related to those of ‘‘sex stuffs’ found in certain algae (Geitler,® 
Moewus’) and Protozoa (Kimball). In Tetraspora lubrica Geitler found 
that centrifugates of culture fluid of one sex cause atypical group-formation 
among cells of the other sex. Copulation, however, does not follow this 
type of group-formation. In Chlamydomonas eugametos Moewus dis- 
covered the same effect of filtrates or centrifugates of culture fluid of one 
sex on individuals of the other sex. In Euplotes patella Kimball found 
that conjugation can be induced between animals of the same mating type 
by placing them in fluid in which another mating type has been living. 

Another possibility is that the phenomena in P. bursaria described in 
this paper are similar to those of the “killers” as reported by Sonneborn!! 
for P. aurelia in view of the following facts: (1) The Ru22 fluid makes the 
animals of the English clone En1 sluggish, darker and distorted in shape and 
these phenomena are similar to those in lethal interactions between diverse 
stocks in P. aurelia. (2) In P. aurelia Sonneborn (unpublished) found that 
one of the ‘‘killers,’’ among other effects, induced pairing as described in 
this paper for P. bursaria. 


* Aided by grants from the Committee for Research in Problems of Sex of the National 
Research Council and from the University of California. 

1 The fluid of two other Russian clones (Ru21, Ru30) appears to have the same 
effects as those of Ru22 fluid. These three clones belong to the same mating type and 
are characterized by a large body and an unusually large micronucleus. These clones 
are vigorous; flourishing cultures can be easily obtained. 

2 Jennings, H. S., and Opitz, P., Genetics, 29, 576-583 (1944). 

3 Jennings, H. S., Ibid., 24, 202-233 (1939). 

4 Of course, it is also possible, though not probable, that the Ru22 fluid causes a 
change in mating type in some of the animals of the English clone, and hence in this case 
it is merely a conjugation between animals of diverse mating types. 

5 In addition to the conjugating pairs, some “‘threes’’ were also observed. (A group 
of three animals in conjugation is to be called a ‘‘three.’’) In each “three,’’ nuclear 
changes occur in all three conjugants and anlagen are formed in the ex-conjugants. 

6 Diller, W. F., Jour. Morph., 59, 11-67 (1936); Caldwell, L., Jour. Exper. Zool., 66, 
371-407 (1933). 

7 As yet only a few clones of these three varieties have been tested. It is likely that 
induced conjugation will be found in more clones when a greater number of clones are 
tested. 

8 Geitler, L., Biol. Zentralbl., 51, 173-187 (1931). 

® Moewus, F., Arch. Protistenk., 80, 469-526 (1933). 

10 Kimball, R. F., Genetics, 27, 269-285 (1942). 
11 Sonneborn, T. M., Amer. Nat., 73, 390-413 (1939). 
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MARKOFF CHAINS WITH REVERSE TRANSITIONS 
By I. OPATOWSKI 
Tus UNIVERSITY OF CHICAGO 
Communicated October 15, 1945 


Consider a system which may take nm + 1 states: 0,1, ..., m. Let 
the probability of transitions (¢ — 1 — 7) and (1 + 1 — 7) during any time 
dt be, respectively, kidt + o(dt) and g,dt + o(dt). Let these be the only 
transitions possible during dt and let the system be in the state 0 at ¢ = 0. 
If Y,(¢) is the probability that the system be in the state 7 at the time ¢, 
then Y,(0) = 1, Y,(0) = O for 2 1 and! 


dY,/dt = RiVi-s — Rieti Yi + 2iVini — gers, (1) 


where by definition Y,4, = Y; = gj; =0ifj7< 0. In certain biophysical 
applications? of (1) the only function of interest is Y,(¢) and k; = k or = 0, 
g: = gor = 0, where k, g are constants > 0. The purpose of this paper is 
to derive, under these assumptions, an explicit expression of Y,(#). The 
method used is that of Laplace transformation which has been introduced 
into this field by H. Bateman.* This transformation seems to be essential 
here because it shows an intimate connection of our problem with Tche- 
bychef polynomials. 

Put yi(s) = fo °e-* Y,(t)dt and consider first the following three processes. 

First Process.—Let ki = RA <t< n), gi: = gOSt<m — 1), kn+1=0. 
Applying the Laplace transformation to (1) we obtain a system of ordinary 


linear equations in y,’s, which gives 


Yn(s) = R"/[En+1 — (R + g)En + kgEn— 1], (2) 
E,, being a continuant, i.e., a determinant |la,, ;||, with 
€1=S+R+ 4, 141 = 2 G-11= hk, a) = 0 for |j — 4 ee 
(3) 
Expanding E,,4; according to the elements of the last colunin we have 
Enti = (Stk + g)En — kgE,-1. (4) 


By (2) and (4) we obtain 
Yn(s) = k"/(sEn). (5) 
Solving (4) as a finite difference equation in E, we have‘ 
Ey = [(kg)"/(o? — DY*[(@ + Va? — 1) — (@ — Vo® = 1)"*1/2, 
where 


o =(s +k + g)/[2(kg)'"] = cos 8. (6) 
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But arccosh ¢ = 16 = log, (o + Vo? — 1) = —log, (o — Vo? — 1). 
Therefore, 


E, = (kg)"/? sin (n@ + 6)/sin 0. (7) 


The factor in @ is here a standard form of the Tchebychef polynomial of 
2nd kind.’ If |c| > 1 we take o = cosh @ and replace in (7) sin by sinh. 
By (7), E, = 0 for 


o = cos (irv) where vy = 1/(n + 1),72 = 1,2, ...,m. (8) 
Therefore, E, considered as a polynomial in s has its m roots s = —k,, where 
k, = k + g — 2(kg)'”* cos (ixv) > 0. (9) 
By (3) the highest term of E, in sis s*. Therefore 
E, = Wri (s + i). (10) 
Consequently by (5) and by known theorems of Laplace transformation :!: * 
Y,,(t)/k" = 1* exp (—kit)* exp (—Ret)*...* exp (—R,t) = (11) 
(Miz B)-! — Dixy (RATGzY jas (By — Be) exp (—Ru), (12) 
Y,(t) = P — (2 Vg/k)!- e-&+0t Y=" He"t/c, (13) 
where * stays for convolution, 1* ... = fo'...dt and 
c, = 1 + (g/k) — 2(g/k)'”* cos (inv), 1/P = Wirt, (14) 
1/H, = 142% ; 2; [cos (éxv) — cos (jrv)], 74 = 2t Whg cos (inv). (15) 
From (14) we have by a known formula’ 
P = [1 — (g/k))/[1 — (g/k)**"). (14’) 


From (13) we see that Y,(t) — P when t > , because (k + g)t — % = 
kt > 0. A simple formula for H; is obtained as follows: from (7) (9) 
(10) and (6) we have 


I=! [cos 6 — cos (jxrv)] = 2- sin (nb + 6)/sin 8, (16) 


which could be derived also from known relations for trigonometric prod- 
ucts.’ Putting 6 = ixv + 6 in (16) we see from (15) that H, is the limit 
for 6 = 0 of 


[cos (ixv + 56) — cos (inv) ]/T14=7 [cos (irv + 5) — cos (jrv)] = 
2[...] sin (trv + 6)/sin [(n + 1)(ixv + 54)], 


where [...] is in both numerators the same. Consequently 


H, = —(—1)!2* sin? (inv). 
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Second Process—Let ky = kR(1<i< n),g,=¢(0<i< n — 2), Rati = 
£n-1 = 0. We have now 


Yn(S) = k"/(sG,) where G, = E, — gEp-1. (17) 
Y,(t) is still given by (11) (12) if the roots of G, are called s = —k, 
(@=1, ..., m). Since? (E,)—0 = k"/P by (10), (9), (14), from (14’) 
we have (G,)<0 = k”. Therefore: 
Rik. . Rn = k*, (18) 
which reduces by one the number of roots to be calculated. We write: 
ki =k +g — 2(kg)”* cos 0, (i = 1, ..., »), (19) 


where the 6,’s are ” values of 6[cf. (6)]. We locate these values so that the 
calculation of the roots may be carried out by routine methods. From 
(7) (17) we have 


' Gy = —(—1)(kg)"/? for 6 = (i — 1)x/n, (i > 2), (20) 
Gn = (—1)*(kg)"/? (g/k)'* for 0 = ix/(n+1),@>1), (20) 
G, = (kg)"/? [nm + 1 — n(g/k)'”] for 0 = 0. (22) 


Therefore, if g/k < (1 + 171)’, so that [...] in (22) is > 0 we have: 
(4 — 1)4/n < 0 < ix/(n + 1), (¢ = 1, ...2), (23) 


which gives a particularly good location of the m roots for large n. If 
[...] in (22) is < 0, we must take (23) for 2 < i < n, which locates only 
n — 1 roots, however, the mth root may be calculated from (18). From 
(18) and (12) we see that Y,(#) ~ last— o. 

Third Process.—Letk; = R(1<i<cn+1),g¢:=2(0<1<"-—1). We 
have now y,(s) = k"/Gn+1. Them + 1 roots of G,+, are located by equa- 
tions (18) to (23) in which, however, m must be changed inton + 1. By 
known methods!: * we have now: 


Y,(t) = k" exp (—Ryt)* exp (—fot)*...* exp (—Ra4il) 
= kort} (ett ji (Rs — B)]-' exp (—Rit). 


Since k, > 0, we have Y,(#) > 0 ast— @. 

More General Processes.—If in the previous processes some g;’s with 
4< n — 2 are zero, it is easy to see’ that the corresponding y, is a product 
of some y,’s of the previous three types, therefore, the Y,,(¢)'is a convolution 
of some Y,’s of the same types. 


1¥For bibliography on the subject cf., e.g., Opatowski, I., these PROCEEDINGS, 28, 
84-88 (1942). The probabilities k,dt+-o(dt) and g,dt+-o(dt) are conditional; they give 
the probabilities of transitions within the interval of time (¢, ¢+d¢) under the assump- 
tion that at the time ¢ the states i—1 and i+1 existed, respectively. 
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2 Opatowski, I., Bull. Math. Biophysics, 7, 161-180 (1945); 8 (1946). 

3 Bateman, H., Proc. Phil. Soc. Cambridge, 15, 423-427 (1910). 

4 Scott, R. F., Messenger of Math., 8, 131-138 (1879). 

5 Szegé, G., Orthogonal Polynomials, New York, 1939, pp. 3, 118. 

6 Cf. Muir, T., Educational Times, 42, 95-96 (1885). 

7 Cauchy, A. L., Algebraische Analyse, Berlin, 1885, pp. 381-382. 

8 Instead of (18) the relation 2} = Tk; =nk+(n—1)g can also be used to calculate one 


of the roots when all the remaining roots are known. 


GENERALIZED ARC-SETS 
By A. D. WALLACE 
DEPARTMENT OF MATHEMATICS, THE UNIVERSITY OF PENNSYLVANIA 
Communicated October 15, 1945 


The purpose of this note is to describe a type of strengthened extension 
set. These sets may be regarded as generalizations of the A-sets introduced 
by W. L. Ayres and G. T. Whyburn and as homotopic analogs of sets stud- 
ied by the latter in his work on cyclic elements of higher order. 

In what follows H will denote a compact (= bi-compact) Hausdorff 
space and subspace will mean closed subset of H. As usual S, is the unit 
sphere of (n + 1)-dimensional euclidean space. A subspace is a C, if all 
mappings of it into S, are inessential (and hence null-homotopic). By a 
T, is meant a subspace all of whose closed subsets are of type C,. Thus, 
in the sense of homotopy, a T, is trivial in dimension 1, at least in so far as 
spheres are concerned. Since an AR can be contracted over itself to a 
point it must bea C, forallm. A p-cellisaT,forp <n. But S, is nota 
C,, and hence not a T,,. 

If M is a subspace then M is termed an extension set in dimension n (or a 
J,) provided that for each subspace X and mapping f: M-X — S, there is 
an f; X — S, such that f = f\M-X. A subset of H is a B, if it is not cut by 
any 7, and is maximal relative to this property. This last concept is due 
to G. T. Whyburn. The sets 7, are homotopic analogs of sets used by 
Whyburn to define cyclic elements of higher order in the same manner in 
which the sets B, are defined. It is to be emphasized that extension sets 
are defined relative to H and the same remark applies to the sets B,. 

The sets J, suffer from the serious defect that they may lie wholly 
within a set of type B,. For example, let N be defined by the equations 
(t>0)x=e'+ |cos t\, y = sin t, let L denote the segment joining (2, 0) 
to (2, —3) and (2, —3) to (0, —3). Finally, if N’ is the image of N in the 
Y-axis and C the circle x? + y? = 1, let H be the union of N, N’, L, L’ and 
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C. Then H is a Bo and Cis a Jo. This latter fact follows readily from 
Theorem 1 below. 

Suppose that X is a subspace and f a mapping of X into S,. Then a 
closed subset Y of H is an essential membrane for f if f can be extended to a 
mapping of any closed proper subset of Y into S, but not to a mapping of Y 
into S,. This concept is due to W. Hurewicz. 

THEOREM I. In order that a subspace M bea J, tt is n.a.s. that for any 
closed set X in M and mapping f : X — Sy, each essential membrane for f be 
contained in M. 

To avoid the difficulty indicated in the example we define a set NV to be 
an A, provided that 

(i) Nisa Jy, 

(ii) If X isa B, and X-N is not a 7, then X is a subset of N. 

THEOREM II. The intersection of any family of A,’s isan A,. Each B,, 
is an Ay. 

This result is well known if m = 0 and His Peanian. Analogous results 
hold for sets involving the higher order cyclic elements. In the proof use 
is made of the fact that J,’s have the same intersection property and that 
each B, is a J,. 

THEOREM III. Let X bea subspace such that H-X is the union of a family 
of pairwise disjoint open sets the boundary of each of whichisaT,. Then X 
ts an A». 

Again when » = 0 and H is Peanian this theorem and its converse are 
well-known results of the cyclic element theory. As to the converse we 
give an example. Let X be the set x? + y? + 2? = 2 and U that part of 
the set (x? + y? = 1,0 < z < 3) which lies outside of X. Then, if H = 
X + U, H — X = U cannot be expressed as the union of a collection of 
pairwise disjoint open sets (other than the family consisting of U itself) 
and U — U is an S; and hence certainly not a 7;. It is easy to see that 
X is a B, and hence a J;. Here H is an ANR and so locally connected in 
all dimensions in the strongest possible manner. 

The space H is said to be of type V provided that for any pair of closed 
sets A and B there is a decomposition of H into closed sets A’ and B’ such 
that A C A’, B CB’, and A’-B’-(A + B) = A-B. If H is metric it is of 
type V but there are spaces of type V which are not metric. The exact 
position of this axiom in the hierarchy of separation postulates is not 
known. It is clearly stronger than normality. 

THEOREM IV. If His of type V then each Ay is an A n+1. 

The proof of this proposition makes use of fact that each 7, isa T,41 and 
each J, is a Jn41. 

It is clear that H is an A, and that any T,, isan A, ifm <n. Also it 
may be shown that if X isa C, and M an A, then M-X isaC,. It follows 
that if His a C, so also is each A, and hence each B,, 
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To generalize the notion of an end-point we say that a subspace E is a 
T,-end-element if each neighborhood of E contains a neighborhood whose 
boundary isa 7,. In this way a point is an end-point if and only if it is a 
To-end-point. 

TuHeEoreEM V. Each T,,-end-element 1s an A». 

By a prime J,(A,) is meant a J,(A,) no closed subset of which is a 
J,(An) without at the same time being a 7,. 

Lemma 1. Let A and B be subspaces such that A,B and A-B are T,’s and 
A+BisaC,. ThnA+ BisaT,,. 

This result is no longer valid if the assumption that A + B be a C, is 
deleted. In general it can only be said that the union of two T;,’s is a 
To+1- 

Lemma 2. If X isan A,,ZaT, and X — Z = U+ V,U|V, then the 
sets U+ Zand V + ZareA,’s. 

With the aid of these results it may be shown that 

THEOREM VI. Jf HisaC,, each prime A, is a B, and conversely. 

THEOREM VII. Let X and Y be subspaces such that X + Y and X-Y 
are A,’s. Then each of X and Y is an Ay. 

For a Peanian H/ the original definition of an A» was as follows: X is an 
A-set if it contains each simple continuous arc whose end-points it contains. 
If now Q,+1 is an (m + 1)-cell it is clear from Theorem I that an A, has the 
property that it contains each Q,+1 whose homology boundary it contains. 
Thus, if in the second example, we add to H the set (x? + y? = 1, 2 = 3) 
then X is no longer an A, since it contains the boundary of a Q, without 
containing the Qo. 

Lemma 3. If HisaT, then H X (01) ts a Ty41. 

With the aid of this result and certain others which follow from it we can 
show that 

THEOREM VIII. Jf HisaC, and not aT, then H contains non-degenerate 
sets B,—1. 

It should be possible to replace the hypothesis that H is a C, by the 
assumption that H/ is C,-1. 

If the mapping f of X into S, is essential but not essential on any closed 
proper subset of X then f is said to be irreducibly essential on X. It may be 
shown, for example, that if H admits an essential mapping onto S, then 
this mapping is irreducibly essential on some subspace. Parenthetically 
we interpolate the following: Suppose that H admits an irreducibly 
essential mapping onto an m-sphere. One might be inclined to conjecture 
that there exists a subspace with the property that every one of its essential 
mappings into S, is irreducibly essential. This need not be true even if H 
is a finite dimensional metric space (of course compact). 

Proposition. Jf H admits an irreducibly essential mapping onto S,, then 
H ts not cut by a T,-;. 
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As of the present writing the validity of this assertion can only be 
affirmed in certain highly restricted situations.} 


1 A detailed discussion of the-sets By, Cy, Jn and T, will appear in the Trans. Am. 
Math. Soc. where pertinent references to the work of Borsuk, Eilenberg, Hurewicz and 
Whyburn will be given. 


TWO MEAN THEOREMS IN HILBERT SPACE 
By Ky Fan 
THE INSTITUTE FOR ADVANCED STUDY 
Communicated November 6, 1945 


The ergodic theory .and the researches concerning the law of large 
numbers in the Theory of Probability have demonstrated the importance 
of the study of the asymptotic behavior of the arithmetical mean 


ty + tet... + Sn 


nN 





of a sequence {x,} of points in Hilbert space. In this note we shall give 
two mean theorems, each of which can be regarded as generalizing the 
following well-known theorem due to J. von Neumann:! For any unitary 
operator U in the Hilbert space and for any given point x in §, the 
arithmetical mean of the n first iterated images of x 


Ux + Ux +... + U% 
n 





converges strongly to a limit point. 
TueoreM I. Let {x,} be a sequence of points in Hilbert space (real or 
complex). If there are two constants ¢, C2 such that the inequalities 


I|xnt1 + Xnp2 + ithe + xn+m||? ane ||. + Xe + eee + xm||? < cm (1) 
and 


|| + xe + eae + Xn+2| |? =“? 2| |x + x2 + mae + xn41|}? 
+ ||or tre +... + x9||? > ce (2) 


hold for all positive integers m, n, then the limit 


lim 21 2+ 1: HX 


n—> ~o nN 





(3) 


exists in the sense of strong convergence. 
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Proof.—We denote as usual by (x, y) the scalar product of two points 
x, yin § and use the following notations: 


So. Se he es A Xa 
((x, y)) = (x, v) + (y, x), 


k 
= a ((%4, Xet+1)), 
n—-1 
Ba = 2 Ox. 
k=1 


Obviously the numbers ((x, y)), ax, 8, are always real. 
It is easy to see that (1) is sae to the inequality 


















































l2 2 Ile 
fs _ sen ff ¢_m_ (lll! salt _ seal, 
—— | - Lt Cy}, 
|| n+ml| n(n-+m) n r+ m n+m T 4 
which can be written 
Is. ae m & 2» ; 
«en x z x,||? 
lace n(n + m) > lel +m ! 
— m B., Bm Bn m 
he — ie +o" Sun) 4) 
m n(n-+m) \n m n+m 
But the hypothesis (1) implies also the boundedness of | [24 P Let ¢3 
be a constant such that 
liga? < @ . 46 = 1,2,3, ...}. (5) 


Then the expression in the first bracket on the right-hand side of (4) is 
less than c, = 2c; + c¢). Hence, we have 














ss _ Sam ||? < i oe (= Bm _ eae) 
||” n + m| n(n + m) m n+m 
or 
Se Setm | _ mM i 2 —_. 
in n + m|\ "Ctaw n+m/ m? 
m Bintimn i Bu 
(n +m)? n (6) 


Now, in virtue of the Schwarz inequality and (5), the numbers a,/k are 
bounded. If we set 


(7) 


= lims 
im sup 7 E’ 
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we have 


i throw. 4 
Bg bkG <a to) (8) 
e(m) being a function of m suitably chosen with lim e(m) = 0. 


Making use of the hypothesis (2) and (5), we can find a positive constant 
Cs such that 


An+1 > ay — cs (k = L, 2. 3 ee a. (9) 
Consequently we have 
m+n—1 ae 
|, pees = Be = be Gy > hy = als 1) cj > N(Qm —_ NCs). (10) 





k=m 


From (6), (8) and (10), it follows that 


2 C4 m r m 2/r 
+ otor es (3+ cn) + () (3 7 cm)) 


a 
(n +m)? 


Sn Snt+m 


Nn n+m 

















(am — nes). (11) 


Similarly, if / is a positive integer not greater than m, we may write the 
inequality” 

















S:__Sutm |? a cae ) 
l n+m ia ty n+m a t + 
seems § i \- se 

( ea 5 + en +m — l) oa (am — Ns). (12) 

According to (11), (12) and 
So sa Sn Sate il’ ols Sata |i" 
=| <Q —- —="} +2) —- "1, 1 
n li iin n+m l n+m (13) 






































we obtain for] S n: 














B-<P+an Gt) <*G Ge) 
in l . 2 ee “idles 7 n+m Rees, 
As 2m 2 2¢4 2(n + m — l) C ) 
e+ mom ncs) + j pater eas rae . + e(l)} + 
Ea Batis 9 . a = 
2(ntn— gt t+ L) er Eger (Gm — Cs). 


(14) 
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If now we keep n, / fixed and let m — o through a suitable sequence of 
integral values such that (a,,/m) — X, (14) becomes 


ae 
< 2ca(1 + 1) + 20m) + 2). (15) 


As this inequality is symmetric with respect to m, /, it remains true even if 
we cease to assume /.< n. This inequality proves the existence of the 
strong limit (3). 

THEOREM II. Let {x,} be a sequence of points in Hilbert space (real or 
complex). If there is a constant c, such that the inequality (1) is fulfilled for 
all positive integers m, n and tf the limit 


fim (eet ae + +. + xntal|? — [lor + me +... + xl? 


n—>@ n 


(16) 





exists, then the strong limit (3) exists also. 

Proof.—Using the same notations as above, we have still the inequality 
(6), which is, in fact, a consequence of the hypothesis (1) alone. But, 
according to the present second hypothesis, we have this time, instead of 
(7) and (8): 


d= lim “ 7 
i k (17) 
and 
B, _ X : 
with lim e(”) = 0. 
By setting 
m Bat m “4 Bm 
et oe (19) 
we infer from (17) that 
lim y(n, m) =X (20) 


for any fixed value of n. 
According to (18), (19), the inequality (6) may be written 


Soto Ce (n)) + 
a ee iy 





Sn Sn+m 





n n+m 











| 
| 





m 27r 
(2 + 4 € “+ cn) — y(n, m). (21) 
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Similarly, for any positive integer / S m, we have 











St Sut m ||? Cs n+m-—l d ) 
i pam? gan er 
sr 2 
eS (+ ein m— 1) — 10, +m — D. (22) 
n+m 2 


Now, using the inequality (13), we can deduce from (21), (22) an in- 
equality which is like (14) and which can be again reduced to (15) by passing 
to the limit m — . Our theorem II is thus proved. 

We terminate this note with the following remark from the point of view 
of the Theory of Probability. If we consider a sequence of chance variables 
instead of a sequence of points in Hilbert space, a similar argument will 
furnish two theorems giving new sufficient conditions for the validity of 
the law of large numbers. The results thus obtained contain as a par- 
ticular case an important theorem due to A. Khintchine,* which asserts 
that any stationary sequence of chance variables obeys the law of large 
numbers. 


1 Neumann, J. von, these PROCEEDINGS, 18, 70-82 (1932). See also Hopf, E., Ergo- 
dentheorie, Berlin, 1937, p. 28. 

2 We use here the inequality anim-1 — lcs = am — mcs, which follows readily from 
(9) and/ S n. 

3 Khintchine, A., Rec. Math. Moscow, 40, 124-128 (1933). 
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